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Abstract 

The authors have employed a numerical procedure to analyze the adhesive contact between a 
soft elastic layer and a rough rigid substrate. The solution of the problem, which belongs to 
the class of the free boundary problems, is obtained by calculating the Green's function which 
links the pressure distribution to the normal displacements at the interface. The problem is then 
formulated in the form of a Fredholm integral equation of the first kind with a logarithmic kernel, 
and the boundaries of the contact area are calculated by requiring that the energy of the system is 
stationary. The methodology has been employed to study the adhesive contact between an elastic 
semi-infinite solid and a randomly rough rigid profile with a self-affine fractal geometry. We show 
that, even in presence of adhesion, the true contact area still linearly depends on the applied load. 
The numerical results are then critically compared with the prediction of an extended version of 
the Persson's contact mechanics theory, able to handle anisotropic surfaces, as ID interfaces. It is 
shown that, for any given load, Persson's theory underestimates the contact area of about 50% in 
comparison with our numerical calculations. We find that this discrepancy is larger than what is 
found for 2D rough surfaces in case of adhesionless contact. We argue that this increased difference 
might be explained, at least partially, by considering that Persson's theory is a mean field theory 
in spirit, so it should work better for 2D rough surfaces rather than for ID rough surfaces. We 
also observe, that the predicted value of separation is in very good agreement with our numerical 
results as well as the exponent of the power spectral density of the contact pressure distribution 
and of the elastic displacement of the solid. Therefore, we conclude that Persson's theory captures 
almost exactly the main qualitative behavior of the rough contact phenomena. 

PACS numbers: 46.55.+d, 68.35.Np, 46.50. +a, 81.40.Pq 
Keywords: contact mechanics, roughness, adhesion, tribology 
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I. INTRODUCTION 



Numerical studies [ij, jsj, j3] have shown that, in case of non-adhesive contacts, when 
an elastic body is brought into contact with a rough surface the true contact area increases 
proportionally to the applied load. To predict such a behavior two main approaches have 
been developed: (i) multiasperity contact theories (originally formulated by Greenwood and 
Williamson (GW) jsj, jof, 0], p], jof) where the contact between the surfaces is modelled 
as an ensemble of randomly distributed Hertzian contacts between the asperities, and (ii) 



Persson's theory of contact mechanics [lOj, [Uj where the probability distribution of the 
contact pressure is shown to be governed by a diffusive process as the magnification at 
which we observe the interface is increased. The scientific community is debating about 
which theory gives the most accurate results. In a previous paper [l2| one of the authors 
(G.C.) has shown that GW-type theories predict linearity only for vanishingly small contact 
areas and load, whereas as the load is increased the theoretical predictions rapidly deviate 
from the asymptotic linearity. This behavior has been shown not to be followed by Persson's 
theory, which predicts linearity between contact area and load up to values of about 15- 
20% of nominal contact area, in agreement with some experimental and numerical results. 
Numerical calculations by Campana et al. [l3(] have shown that Hertzian-type regime, which 



is the basis on which GW and similar theories have been developed, occurs only at relatively 
small loads, thus indicating the inadequacy of GW-type theories at higher loads. 

As already observed the original version of Persson's theory describes the interfacial con- 
tact pressure through a parabolic partial differential equation where the diffusivity term is 
calculated under the approximation that the Power Spectral Density (PSD) of the elastically 
deformed surfaces is equal to the PSD of the underlying rough surfaces [l^, [l5[, i.e. assum- 
ing that the diffusive term that one would obtain in case of full-contact conditions remains 
exactly the same also in case of partial contact conditions. The stored elastic energy (or, in 
case of sliding of viscoelastic solids, the friction coefficient) is, instead, calculated assuming 
that the PSD of the deformed surface is the product of the PSD of the underlying rough 
surface times the fraction of contact area at the given resolution. This, in particular, can be 



shown to be coherently derived by the theory itself, see Ref. [15|. Of course in full contact 
conditions Persson's theory is exact, but in case of partial contact it has to be possibly ver- 
ified. For this reason, there is not yet a clear evidence about the correctness of the factor of 
proportionality between contact area and load predicted by Persson's theory. Indeed, there 



are numerical investigations [3| , [IJ] of non-adhesive contacts between rough surfaces, which 



show that Persson's theory underestimates the contact area, although its main qualitative 
prediction seems in very good agreement with numerical calculations, as proved in Ref. jl3| . 
where Green Function Molecular Dynamics (GFMD) numerical calculations have been em- 
ployed to show that the PSD of the deformed surface has the same power law exponent as 
predicted by Persson's theory [lEj 

In this paper the authors make an attempt to give an additional contribution in this 
direction. We extend the analysis to include adhesive interactions at the interface of the 
contacting bodies, which become more and more important as the length scale of observation 
is decreased and may dominate the contact behavior of micro- and nano-mechanical and 
biomechanical systems. We indeed focus on the adhesive contact between a semi-infinite 
half space and a randomly rough surface with roughness in only one direction, and compare 
our numerical predictions with the results of Persson's theory. We observe that the original 



version of Persson's theory [10[, [Uj was conceived to deal with isotropic surfaces, but in our 
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case the surfaces is strongly anisotropic, being rough only in one direction. Therefore, in 
order to compare theoretical and numerically calculated data we have employed an extended 
version of Persson's theory |16|], which is able to handle anisotropic surfaces as in the case of 
ID rough surfaces. There are mainly two reasons for studying a ID rough surface: (i) first of 
all one should consider that surface roughness is characterized by a large number of length 
scales, which can cover 3-4 decades and even more. Therefore, in order to get physically 
meaningful results, one needs to include all the spectral components of the surface roughness 
in the analysis. However, increasing the number of length scales rapidly increases the number 
of points where the numerical solution has to be sought, and, in turn, the computation time. 
This problem is strongly reduced in case of ID roughness so that one can include in the 
analysis more than 3 decades of length scales; (ii) secondly we must also observe that rough 
surfaces, encountered in many practical applications, are often strongly anisotropic mainly 
as a result of machining and surface treatments (e.g. unidirectional polished surface which 
present wear tracks along the polishing direction, although the resulting roughness is not 
strictly ID). Thus, from a practical point of view, it is also very important to test Persson's 
theoretical prediction for anisotropic surfaces. 

In the last years scientists have been developing ad hoc numerical methods to treat the 
problem of contact mechanics between randomly rough surfaces. Here we would like to 
recall the methodology proposed by Robbins and co-workers [3|, who developed a Coarse- 
Graining FEM (CGFEM) approach, and that conceived by Campafia and Miiser [l^], who 
have developed a Green's Function Molecular Dynamics (GFMD) approach to deal with 
such a problem. Here we employ a different methodology to deal with adhesive contact. 
The methodology, already presented by one of us in Ref . [iSj , is based on a pure continuum 
mechanics approach and belongs to the class of Boundary Element Methods (BEM), since 
it also makes use of Green's function to solve the problem. This allows us to reduced the 
problem to Fredholm integral equation of the first kind with a logarithmic kernel. We stress 
that the position of the edges of each contact patch is not known a priori and must be 
determined by requiring that the total energy of the system is stationary, i.e. we are dealing 
with a free-boundary problem. 

The numerical procedure has been designed in such a way to never loose resolution 
even when the single contact spot is below the smallest length scale of observation. We 
show also that the numerical complexity of the problem can be strongly reduced since the 
thermodynamic state of the system only depends on the size of each contact area and on 
the pressure distribution in the contact area. Thus, in our boundary element approach, 
only the contact patches need to be discretized (we use an ad hoc adaptive grid) and the 
solution of the Fredholm equation, which is obtained by means of matrix inversion, has to 
be determined only for a limited number of points, i.e. only those belonging to the true 
contact area. 

II. THE NUMERICAL MODEL 

We consider a periodic contact as shown in Figure 1 where an elastic layer of thickness d is 
interposed between a flat rigid plate (upper surface) and a periodically rough rigid substrate 
with wavelength A (bottom surface). We assume that the rough surface has roughness in 
only one direction and is smooth in the orthogonal direction. Under these conditions the 
problem at hand is a periodic plane problem, i.e. the stress, displacement and strain fields 
only depends on the x and y coordinates shown in Fig. 2 and are periodic functions of 
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FIG. 1: An elastic layer of thickness d in adhesive contact with a rough periodic substrate of 
wavelength A. 
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FIG. 2: The definition of substrate displacement Utot; elastic layer average displacement u„ 
substrate penetration A. 



and 



period A. Fig. 2 shows, in particular, the total displacement Utot of the substrate, the 
average displacement Mm of boundary of the deformed layer and the penetration A of the 
rigid substrate into the elastic slab. These three quantities are shown to satisfy the following 
relation 

Mtot = A + Mm (1) 

We will focus on the pressure distribution a (x) and the displacement u {x) of the elastic 
solid at the interface. In Refs. [3] and [l^ G.C has shown that the unknown pressure 
distribution in the contact area fl can be determined by solving the following Fredholm 
integral equation of the first kind with a logarithmic kernel as 



[X 



s) a (s) ds = [h (x) — /^max] + X E Q 



(2) 
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where fl = [a,, 6j] is the unknown contact domain to be determined as shown below. 
The quantities a, and bi are the unknown coordinates of i-th contact patch with < bi and 
i = 1,2, L, where L is the unknown number of contacts. In Eq. ([2]), assuming the elastic 
layer is infinitely thick (i.e. d — +cxd), the kernel is 
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and represents the Green's function of the semi-infinite elastic body under a periodic loading, 
i.e. it represents the displacement u (x) — Mm caused by the application of a Dirac comb with 
peaks 5 {x — n\) separated by a distance A. Here E and z/ are the Young modulus and the 
Poisson ratio of the elastic layer. In Eq. ([2]) and the quantity h {x) represents the heights of 
the rough profile measured from its mean plane. Since we are considering a periodic problem 
h ix) can be written as Fourier series 



-oo 



h {x) = 2_^hm COS {mqox + 0^) (4) 

m=l 

where the fundamental wave vector is go = 27r/A. Also we have defined in Eq. ([2]) the 
quantity /imax = max [h (x)], which is the maximum height of the substrate roughness. Once 
the pressure distribution is known the elastic displacements at the interface can be easily 
determined through the equations 

u{x) — u^ = — I Q {x — s) a {s) ds; x G D — ^7 (5) 

u{x) — Uja = h (x) — /imax + ^5 X E Q 

where D = [—A/2, A/2]. Of course for a infinitely thick layer {d — )■ +00 as in our case) the 
average displacement is also infinitely large except when the u = 0.5, but the differ- 
ence u (x) — Um is always finite [3], [l9|, and can be interpreted as the additional elastic 
displacement of the solid due to the presence of roughness at the interfaces. In order to 
close the system of equations we need ad additional condition to determine the yet unknown 



contact domain Q. To this end (see also Ref. 18j), we first observe that for any penetration 



A, we can calculate the pressure distribution at the interfaces through Eq. ([2]), and the 
interfacial elastic displacement through Eq. as functions of the unknown coordinates 
Cj and bi of the i-th contact area. To calculate the exact values of the quantities and 
bi, given isothermal conditions, we need to find the stationary point of the free interfacial 
energy f/tot (oi, &i, a^, bi. A) of the system for a fixed value of the penetration A, this is 
the same as requiring that 



dtti / A ), \ dbi 

The interfacial energy (see Ref. [l8|) is 

t/tot = t/el + t/ad (7) 
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where we have defined the interfacial elastic energy Ud as the amount of elastic energy stored 
in the solid as a consequence of the elastic deformations caused by the substrate asperities 
i.e. 

Uei{ai,bi,...,aL,bL,^) = -^ a (x) [h (x) - hrai,^ + A] dx. (8) 

^ i=i 

and the adhesion energy is 

U^d{ai,bi,---,aL,bL) = \l ^ + [h' {x)fdx. (9) 

where 7 is the Dupre energy of adhesion per unit area. Eqs. ([2]), ([5]) and ([6]) constitute a set 
of closed equations which allows, for any given penetration A, to determine the coordinates 

and bi of each contact patch, the pressure distribution at the interface, and all other 
thermodynamical quantities. For the numerical implementation the reader is referred to 
Ref. [l8|, here we just describe some numerical techniques which are peculiar to the problem 
we discuss in this paper. Let us assume that we know the solution of Eqs. (E]) for a given 
penetration A: the knowledge of the contact region Vt = [jf^i[ai,bi] is sufficient to fully 
characterize the system. Eq. ([2]) determines the stress we solve it iteratively, through 

the Gauss-Seidel algorithm. Special care has been taken to guarantee stability and accuracy 
of the solution by choosing a suitable sampling grid: the domain Q is discretized in a 
non-uniform, adaptive way, so to ensure that more points are employed close to the edges 
of the contacts, where the stress distribution, because of adhesion, presents a square root 
singularity. As long as the stress is known, the deformed profile of the elastic layer follows 
from Eq. ([5]). The interfacial energy is then given by Eqs. ([7]), ([8]) and IQ. 

In summary, at a low level in our numerical implementation of the algorithm there is 
the Solver, a software code that, given the penetration and the contact regions, calculates 
everything else. On top of it we built another piece of software in charge to adjust the 
position of the contacts boundaries ai, 61, . . . , a^v, ^at so to minimize the interfacial energy. 
We employed a conjugate gradient method in the version given by Polak and Ribiere [20|. 
Unfortunately the problem is more complex than a minimization in a 2A^-dimensional space. 
Starting form an arbitrary configuration, some of the contact boundaries can acquire the 
same value, meaning that either a contact is detaching or two contacts are coalescing together 
into a single bigger one. If this happens, the minimization has to be restarted in a different 
number of dimensions. 

Furthermore, a constraint must be accounted: in principle we can determine the con- 
figuration corresponding to any penetration and contacts, but we must impose that in the 
non-contact regions the elastic layer never intersects the substrate. For instance, if we con- 
sider a physical configuration minimizing the interfacial energy, and then we increase the 
penetration pushing the substrate against the elastic layer, the starting point for the new 
conjugate gradient minimization may show an intersection between the elastic layer and 
some peaks of the substrate in the non-contact regions. This indicates that a new contact 
has to be added before starting the minimization. A specific procedure inside our software 
is in charge to detect all the intersections between elastic layer and substrate, so to enforce 
the physical constraints of the problem. Given the penetration A, the search of the contact 
domain Q that minimizes the interfacial energy is challenging: not only the 2A^ variables 
ai, 61, . . . ajv, &Ar are unknown, but also the number N of contact regions is unknown! The 
solution of the problem resorts to conjugate gradient minimization alternated to searches 
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for intersections between the elastic layer and the block. The minimization procedure stops 
when the conjugate gradient ends successfully, i.e. it is not interrupted by a coalescence of 
detachment of contacts, and the successive search for intersections confirms that there are 
no intersections in the non-contact regions. Although we took special care to guarantee that 
the minimization procedure would converge for moderately large variations of penetration, 
we observed that the most reliable approach involves many small increments of penetration 
starting from (non-contact) up to desired value, while optimizing of the solution at every 
intermediate step. The solution of the contact problem in presence of adhesion is not unique, 
that is, for the same penetration more configurations are possible depending on the loading 
history. Nonetheless the contact pattern occurring with increasing penetration is uniquely 
identified and it is the most suitable solution to represent the non-adhesive contact in the 
limit of vanishingly small adhesive bonds. As a final remark, we observe that this algorithm 
cannot be used to solve the problem without adhesion: in this case the solution is still a 
stationary point satisfying Eq. ([6]), but the profile of the elastic layer in the non-contact 
regions is always tangent to the substrate near the crack tips and bi. An infinitesimal 
motion of any of the crack points decreasing the contact region would cause an intersection 
between elastic layer and substrate. In other terms, the solution lies always on the boundary 
of the subset of identified by the physical constraint of non-intersection. The solution 
would not be a minimum without such a constraint. 



III. PERSSON'S THEORY FOR ANISOTROPIC SURFACES 

The aim of this paper is to compare the numerical results with analytical predictions 
of one of the most promising and also strongly debated theory of contact mechanics, i.e. 



the recent theory by Persson [10|, However, the surface we are considering is strongly 
anisotropic, in fact it is rough in only one direction and smooth in the orthogonal direction. 
The original theory proposed by Persson was, instead, conceived and developed for perfectly 
isotropic surfaces. Therefore in order to carry out the analysis we need to extend this theory 
to the case of anisotropic surfaces. This extension has been obtained in Ref. [l6|. Here we 
briefly summarize the main equations. Persson's theory removes the assumption, which is 
implicit in the multiasperity contact theories, that the area of real contact is small compared 
to the nominal contact area. On the contrary, Persson focuses on the probability distribution 
P {a, () of normal stresses at the interface, which depends on the magnification at which the 
contact interface is observed. To calculate the governing equation of P {a, (), Persson moves 
from the limiting case of full contact conditions between a rigid rough surface and an initially 
fiat elastic half-space [lo[. In such conditions the PSD of the deformed elastic surface is equal 
to C2D (q) where C2D (q) = (27r)"^ / (fx {h (0) h (x)) e-^^'^ is the PSD of the rigid rough 
substrate (the quantity h (x) is the rough substrate height distribution, x = (x, y) is the in- 
plane position vector and the symbol () stands for the ensemble average). Considering that 
for a perfectly elastic material the elastic modulus is frequency independent, it can be easily 



shown (see Ref. [16|]) that even for the general case of anisotropic surfaces Persson's theory 



states that the stress probability distribution P (cr, () must satisfy the following relation. 

where the magnification C, = q/qo, and a is the interfacial stress in the apparent contact 
area at the magnification (. The diffusivity function / {() = q^G' (q) cTq, where (Tq is the the 
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average normal stress in the contact area and G (q) is calculated in full contact conditions 



as 

1 / E V 1 



G'(?) = ^(rf^) ^([VMx)f), (11) 



'0 



where 



{[Vh{^)]\ = I d\' q''C{c{) (12) 
-'|q'l<g 

is the mean square value of the slope when the surface is observed at the magnification 
C = q/Qoj that is when all harmonic components of the spectrum with wavevector above q 
are filtered out. Then, Persson assumes that Eqs. ffTUl) and ffTTl) hold true also in partial 
contact (this is of course an approximation since the PSD of the deformed surface in partial 
contact conditions cannot be the same as that of the rigid rough substrate) and to account 
for partial contact the following initial and boundary conditions are enforced 

P{a,l) = 6{a-ao) 

P{-cT,, = (13) 
P(oo,C) = 

Here cxa (C) is the tensile stress needed to cause detachment over a strip of length 27r/g, 
recalling the Griffith criterion in plain strain 22| we get 



^a(C) 



^2 



E*ics (C)g 



.1/2 



(14) 



where E* = E / (1 — t>^). In Eq. (1131) the quantity 7efr (C) is the apparent energy of adhesion 
at the interface defined as [lH, |2l| 

- 7eff (C) A (C) = t/ad (C) + (C) (15) 

where A {Q is the apparent contact area at the magnification C,. The calculation of the 
elastic energy Ue\ (C) must take into account that because of partial contact conditions the 
elastic energy stored at the interface is less than what would be stored in case of full-contact 
conditions. This is necessarily true because only where contact occurs the elastic solid 
conforms the underlying substrate, whereas outside of the contact regions the elastic surface 
is much less deformed. Persson accounts for this fact by assuming that, when calculating 
the interfacial elastic energy Ueu the PSD of the deformed (initially flat) elastic surface is 
equal to C2D (q) A (g) /Aq where A (g) /Aq is the fraction of apparent contact area at the 
length scale 27r/g i.e. 

f/ei (C) = -E* [ d'qqC2B (q) A (q) (16) 

This result, that at the beginning was only conjectured by Persson [lo|, recently has been 
demonstrated to directly follow from Persson's theory itself [3]. Analogously the adhesion 
energy Uad (C) is calculated as 

POO 

f/ad (C) = -lA (Cm..) / dx{l+ ex) e-^ (17) 

^0 
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where 

e = [ d'qq'C,^ (q) 

|q|>C<?o 

and Cmax = Qi/Qo, Qo = 2tt/X, and qi = 2tt/\i where Ai is the shortest length scale of the 
rough surfaces. Equations (fTOl [T3l [HI [151 [ID and [17]) can be solved to calculate the stress 
probability distribution P (a, C) and hence the apparent contact area A {() /Aq as a function 
of the magnification ( [ll|, [2l| 



MO 

A. 



/oo 
P{aX)da (18) 



The separation s = /imax — A can be calculated observing that the change of total interfacial 
energy Utot = Uei + U^a must be equal to the work done by the applied load, i.e. 

dUtot = (ToAodA = -(ToAods (19) 

which gives 

s= / /yr ^cr- (20) 

In case of non-adhesive contact i.e. when 7 = 0, Persson has shown that 

and Eq. ( ITS]) simply becomes 

A{0 



.0 -^VG(C)j ^''^ 

The above formulation holds true also for anisotropic surfaces, in particular if the sub- 
strate has roughness in only one direction, e.g. along the x-direction, we get ([V/i(x)]^)^ = 

(^{dh/dx)'^) ^. In such a case we have h{x.) = h{x) and one can easily show that 

C2D(q) = C{q,)6iqy) (23) 

where C (g) = {2n) ^ J dx (h (0) h (x)) e"*"^^ is the PSD of the x-profile of the surface. Using 
Eq. ([23]) one simply obtains 

{{dh/dxf)^ = r dq,qlC{q,) (24) 

J -q 

Eq. (|T8]) gives the apparent contact area at the resolution A (g) = 27r/g as a function of the 
applied load (Tq. However, we are interested in calculating the real contact area, which can 
be obtained by replacing C, with Cmax- 
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IV. ROUGH PROFILE GENERATION 



In order to carry out the numerical simulations and compare the results with the theoret- 
ical predictions, we need to numerically generate a rough profile. We have opted for a fractal 
self affine rough profile. For any self affine fractal profile h (x) the statistical properties are 
invariant under the transformation 

x^tx; h-> t^h (25) 
in such a case it can be shown that the PSD of the profile is 

C (q) = Co i^-^j (26) 

where H is the Hurst exponent of the randomly rough profile, and is related to the fractal 
dimension Df = 2 — H. In order to carry out the numerical calculations we have utilized 
a periodic profile with Fourier components up to the value qi = Nq^ (in this case Ci = ^)- 
However we need to determine the amplitudes hm and the phases (pm of the harmonic terms 
[see Eq. (jlj)]. It can be shown that in order to satisfy the translational invariance of the 
profile's statistical properties (which implies that the autocorrelation function satisfies the 
relation {h (x') h{x' + x)) = {h (0) h{x))), it is enough to assume that the random phases 
(pm are uniformly distributed on the interval [— 7r,7r[. In such a case the autocorrelation of 
the profile takes the form 



{h {x') h{x' + x)) = COS {mqox) (27) 



m=l 



Now we need to calculate the quantities (hf^). To this purpose let us calculate the PSD of 
the periodic profile of Eq. (jl]) . By using the definition we get 

^ /• /A2 \ N 

C (q) = {2nr' I (^^0^) = Ei K^^> ^ - ^^o) + (hi) 5 (q + mqo)] 

m=l m=l 

(28) 

from which it follows that 

C {-mqo) = C {mqo) = C^ = ^6 (0) 

Using Eq. (12^ and observing that Co = (hf) 6 (0) /4, one then obtains 

K) = {hl)m-(^^-^^^ (29) 

Hence, the quantity {h"^) can be determined once known {hf) and the Hurst exponent of 
the surface. Now observe from Eq. (127|) that (^h (x)^) = J2m=i (^m) Z^; and using Eq. ( 129|1 
one obtains 

<M^)^>=^i:--^^^^^^ (30) 



•m=l 
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FIG. 3: The deformed shape of the elastic body at three different separations, s = 86/im (blue), 
s = 48/um (red), and s = 9.6fim (green). The rough rigid substrate profile is shown in black color. 

Therefore, if one knows the rms roughness of the surfaces /imis = can calculate 

(hi) and therefore all the other quantities {h"^). However to completely characterize the 
rough profile we still need the probability distribution of the amplitudes h^- There are 
several choices, however the simplest assumption, as suggested by Persson et al. in Ref. 
(23| . is that the probability density function of is just a Dirac's delta function centered 
at [4CJ6{0)f^ ^ 2V2^CU^/L^/^, i.e. 



p{h„,) = S{hm-2J-C'J'\ (31) 




where we have used that (5 (g = 0) ~ L/ (27r). In can be shown [23] that this choice guarantees 
also that the random profile h (x) has a Gaussian random distribution. 



V. RESULTS 

We assume that the elastic block is a soft perfectly elastic material with elastic modulus 
E = IMPa and Poisson's ratio u = 0.5, i.e. we assume that the material is incompressible. 
We assume also that the change of surface energy upon contact between the two surfaces 
(i.e. the Dupre energy of adhesion) is 7 = 0.03J/m^. Calculations have been carried 
out for 11 different realizations of a rough self-affine fractal ID profile. The profile has 
a fractal dimension Dj = 1.3 (i.e. the Hurst coefficient is H = 0.7), with root mean 

square roughness (/i^)^^^ = 50yum. The self-affine profiles have spectral components in the 
range qo < q < qi. We have used A = li^jq^ = 0.01m and qi = lO^go- The numerical 
calculations have been carried out for different values of the separations s = /imax — A, 
which is defined as the distance between the mean plane of the deformed surface and the 
mean plane of the rough surfaces. In Fig. [3] we show three different shapes of the deformed 
profiles at three different values of the separation: s = 86/im (blue), s = 38/xm (red), and 
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FIG. 4: The PSD C„ (q) of the rigid substrate profile (black solid line) compared to that of the 
deformed shape of the elastic body for two different dimensionless loads ctq/E* = 0.005 (blue), 
and (Tq/E* = 0.020 (red). Solid lines refer to numerical calculations whereas dashed lines refer to 
Persson's theory. The agreement between Persson's theory and numerical calculated predictions 
is qualitatively very good. Indeed the slope of the PSD predicted by Persson (in the mid range 
of g-vectors where the influence of adhesion is negligible) is almost the same as the numerically 
calculated one (see text): in both cases the PSD of the elastically deformed solid is Cu (q) ~ 



s = 9.6/im (green). The black line instead represents the rigid rough substrate profile. A 
deeper analysis of the figure shows that, not depending on the separation s, full contact 
occurs between the elastic block and the short wave length corrugation of the rough rigid 



profile. This is in agreement with some theoretical arguments [Uj, [21[ which predict 
that this situation should occur when the Hurst exponent of the rough profile is larger 
than 0.5, as in our case. This is also confirmed in Fig. H] which shows the PSD of the 
rough surface and that of the deformed elastic surface as a function of the wave-vector 
q (in a log-log diagram) for two values of the applied stress ao/E* = 0.005 (blue), and 
(Tq/E* = 0.020 (red), where E* = E / {1 — z/^). Indeed, we observe that for large g- vectors the 
PSD Cu{q) = {2tt)~^ J (ix (m (0) M (x)) e"*"^"* of the numerically calculated deformed profile 
u (x) becomes almost perfectly parallel to the PSD C (q) of the rigid substrate: this means 
that the spectral content of the deformed body profile at short wavelength is just the same as 
that of the rough rigid profile, and therefore that full contact occurs between the elastic body 
and the substrate at short wavelengths. Fig. HI shows also, as expected, that, as the load 
is increased, the quantity Cu (q) continuously approaches the PSD C (g) of the rigid rough 
profile and must become equal to C (q) at very high loads, i.e. when full contact occurs. In 
Fig. Hlwe also compare the numerically calculated PSD of the deformed surface (solid line) 
with Persson's theoretical predictions (dashed line). We first observe that there is a non 
negligible shift between Persson's results and our numerically calculated ones. However the 
two curves run almost perfectly parallel especially in the mid range of wavevectors where 
the best fit of our numerical results gives C„ (g) ~ g"^ ''^. The value —2.73 is very close to 
the value —2.7 predicted by Persson. Indeed, Persson's theory relies on the assumption that 
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FIG. 5: The quantity (q) /E*'^ for full contact conditions (black solid line) and for two different 
dimensionless loads ctq/E* = 0.005 (blue), and ctq/E* = 0.020 (red). Solid lines refer to numerical 
calculations whereas dashed lines refer to Persson's theory. Also in this case the agreement between 
Persson's theory and numerical calculated predictions is qualitatively very good, at least in the 
mid range of (^-vectors where the influence of adhesion is negligible: both calculations predict 
a (q) « q-^. 




FIG. 6: The true contact area A (Cmax) in units of the the nominal contact area as a function of 
the dimensionless applied load ctq/ E* . Numerical predictions are in red, whereas Persson's theoret- 
ical calculations are in black colour. Observe that both numerical calculations and Persson's theory 
predict linearity between contact area and load. However Persson's theory predicts a coefficient of 
proportionality which is about half of the numerically calculated one. The dashed line represents 
Persson's calculations in absence of adhesion interactions. Notice that because of the low amount 
of adhesion the curve does not differ significantly from the solid line. 
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FIG. 7: The separation s in units of roughness /irms as a function of the dimensionless applied 
load (Tq/E* in a log- linear scale. The agreement between numerical data (black curve) and the 
Persson's theoretical predictions (red curve) is almost perfect. The numerical calculations deviate 
from theoretical predictions only at small loads, since our system is a finite system so a finite 
value of the separation s necessarily exists at which contact area goes to zero and therefore also 
the applied load. Persson's theory instead has been developed for infinite systems. In this case 
indeed the rough surface has arbitrarily many and arbitrarily high asperities, which always allows 
the contact between the two solids to occur for arbitrarily large surface separations. 

the PSD of the deformed surfaces Cu (q) = C (q) A (g) /Aq, where C (g) is the PSD of the 
rough substrate. Now, in case of self affine fractal surfaces, using Eq. ( |26|) and observing 
that if one neglects adhesion (this is correct in the mid range of wavevectors q) Eqs. (ITTjl 
and fl22|) give A(q) /Aq ^ q^~^, one obtains from Persson's theory that Cu{q) ~ q^^'^^^\ 
Being, in our case, H = 0.7, we have C„ (g) ~ g^^'^in perfect agreement with our numerical 
calculations. Same conclusions can be found if one observes Fig. E] which shows in a 
log-log diagram the power spectral density Co- (g) = (27r)~^J dx {a (0) a (x)) e~^'^^ of the 
stress distribution at the interface in units of E*'^. Of course this is not unexpected since the 
Ccr (g) and Cu (g) are related each-other through (g) = \E*'^q^Cu (g) [l5[, [l6[. Using that 
Cu (q) ~ g~'-^^^\ and assuming the adhesion interaction is not important (which occurs in 
the mid range of wavevectors g), one obtains, (g) ~ q~^ . This is indeed confirmed in 
Fig. [5] which shows that Persson's prediction (dashed lines) and numerical calculations (solid 
lines) run parallel to each other in low-mid range of g-vectors. However, as g is increased 
the numerical calculated PSD (g) rapidly changes its slope. This, in turn, becomes 
almost equal to that of the PSD of interfacial normal stress distribution that would be 
obtained in full contact conditions (solid black line in Fig. [5]), and confirms, what we have 
already observed above, i.e. that because of adhesion the short wavelength roughness of 
the underlying rigid profile is in full contact with the elastic body. Figure [6] shows the true 
contact area A ((^max) Mo vs. the dimensionless load ctq/E* calculated through Persson's 
theory for adhesive contact and the one computed by our numerical code. Figure [6] confirms 
the linearity between contact area and load, however it also shows a significant disagreement 
between Persson's theory and our numerical calculations. In particular, Persson's theory 
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FIG. 8: The probability function -P (c, Cmax) of interfacial pressure distribution a{x). Points are 
numerical predictions whereas dashed lines are Persson's results. The trend is qualitatively the 
same, although it is quantitatively different. The reason for such a difference is that numerical 
calculations and Persson's theory do not predict the same value of the contact area for any given 
applied load (see Fig. [6]). Notice that, the tail of numerically calculated probability distribution at 
negative loads is an effect of the adhesion interaction which has been introduced only through the 
surface energy, i.e. with the inclusion of an interaction force with an infinitesimally short range. 
We also present a best fit based on double Gaussian probability distribution (see text). 



predicts a contact area which is about 50% less than that calculated with our numerical 



code, in agreement with some molecular dynamics simulations [Ij], where Persson himself 
has found a difference between the numerically calculated contact area and the theoretical 
value of about 30% for a two-dimensional rough surface. The dashed curve in Fig. E] 
represents Persson's theoretical predictions when adhesion is not included in calculations. As 
expected adhesive interactions lead to an increase of the contact area. However, being in our 
case the amount of adhesion energy relatively small, this effect is only marginal. We observe 
that the large discrepancy between Persson's predictions and our numerical calculations may 
be explained, at least partially, by the fact that Persson's theory is a mean field theory in 
spirit, so it should work better for 2D rough surfaces rather than for ID rough surfaces. 
Fig. [7] shows the dimensionless separation s / hj-^^ as a function of the dimensionless load 
ao/E*. Numerical predictions (black line) are compared to Persson's ones. The agreement 
is very good except at lower applied loads. Indeed at small loads the numerically calculated 
separation drops off faster than predicted by Persson's theory. The same effect has also been 
observed in molecular dynamics calculations The explanation for this behavior is that 
numerical calculations have been carried out for a finite system, whereas Persson's theory 
is for an infinite system. Indeed an infinite system has many (arbitrarily) high asperities, 
which always allows the contact between the two solids to occur for arbitrarily large surface 
separations. But a finite system has asperities with height below some finite length /imax, 
and for u > /imax no contact occurs between the solids and ctq — 0. Fig. |8] shows the 
probability density function P {a, Cmax) of interfacial normal stress distribution a [x) at the 
highest magnification. We observe that Persson's predictions and numerical data agree only 



15 



OF 




— 10 i J I I i i i I i i i I i i i ^_ 

-0.2 0.0 0.2 0.4 0.6 

0-/E* 



FIG. 9: The logarithm of the probability function -P (o", Cmax)- Points are numerical predictions 
whereas dashed lines are Persson's results. We observe that the tail of the probability distribution 
at large values of a follow exactly a Gaussian distribution, whereas the tail obtained for negative 
value of a is not Gaussian. 

at a qualitative level, but strongly deviates from a quantitative point of view. The reason of 
this quantitative disagreement can be easily understood if one recalls Eq. f llSp which states 
that the integral of the P {a, Cmax) must be proportional to the true contact area A (Cmax) 
and considers that Persson's theory predicts a contact area smaller by a factor 1/2, if 
compared to numerical calculations. We also present with solid lines the best fit obtained 
with the double Gaussian distribution 
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P (a, Cmax) = —m e-(--'^"+-)'/^« - e-('^+-o+-)V4G (32) 

2 (ttG) 

where we have relaxed the quantities cXa and G. Fig. |H] and even more M show that, at 
least when the amount of adhesion energy is small (as in our case), Eq. f l5^ is a good 
approximation of the numerically calculated stress probability function P (a, Cmax)- Indeed 
Fig. |9] shows that, at high values of a, the tail of the stress probability density function 
is Gaussian. Notice that, the tail of numerically calculated stress probability distribution 
at negative loads is an effect of the adhesion interaction which has been introduced only 
through the surface energy, namely by means of an interaction force with an infinitesimally 
short range. This even allows that infinite negative values of a can occur at the interface. 
We also observe that the value of calculated by fitting the interfacial stress probability 
density function differs from that calculated in the spirit of Persson's theory through Eq. 
dH]). However, it is possible to show that a better estimation of the tensile stress cTa, can 
be obtained by assuming in Eq. f|T^ that the width of the detached region is a factor 1/2 
smaller than that originally assumed in Persson's theory fv^. In such a case the estimated 
tensile stress would increase of roughly a factor 1.4 thus making cxa closer to the value we 
have found in Fig. 13 
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VI. CONCLUSIONS 



The authors have carried out detailed numerical calculations to determine the contact 
area, the stress distribution, the penetration and elastic deformation of an infinitely thick 
layer in adhesive contact with a rough strongly anisotropic rigid surface. The numerical 
predictions have been compared in detail with those of an extended version of Persson's 
theory able to deal with adhesive contact between anisotropic rough surfaces. It is shown 
that, for any given load, the value of true contact area predicted by Persson's theory sig- 
nificantly differs from the numerically calculated ones, the first being smaller by a factor 
PS 1/2. This may also depend on the fact the Persson's theory is of the mean-field type and, 
therefore, should work well in higher dimensions than in ID. However, the predicted value of 
separation matches almost perfectly the numerical data. We have also compared the power 
spectral density of the deformed elastic surface and of the stress distribution at the interface 
as obtained by numerical calculations and Persson's theory. We observe that both theory 
and numerical calculations predict the PSDs to follow a power law with almost the same 
exponents. This extends to the case of adhesive contact, what has been found previously 
for adhesionless contact by other authors. However, we also observe, in agreement with 
Persson's theory of adhesive contact, that at high magnification the exponent of the power 
law changes in such a way to suggest that the elastic solid lies in full contact condition 
with the fine microstructure of the rough surfaces. We conclude that, Persson's theory is 
able to capture the main physics behind contact mechanics of rough surfaces independently 
of whether adhesion is present or not. However, we also observe that from an engineering 
point of view a better estimation of the contact area would be very useful in practical ap- 
plications as in case tires, mixed lubricated interfaces, and seals, where the kinetic friction 
or the amount of leakage should be accurately predicted for design purposes. Therefore an 
improvement of the theory, which allows to better estimate the real contact area, would be 
strongly appreciated by engineering community. 
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